Quantum anisotropic Heisenberg chains with superlattice structure: a DMRG study. 
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Using the density matrix renormalization group technique, we study spin superlattices composed 
of a repeated pattern of two spin- 1/2 XXZ chains with different anisotropy parameters. The mag- 
netization curve can exhibit two plateaus, a non trivial plateau with the magnetization value given 
by the relative sizes of the sub-chains and another trivial plateau with zero magnetization. We 
find good agreement of the value and the width of the plateaus with the analytical results obtained 
previously. In the gapless regions away from the plateaus, we compare the finite-size spin gap with 
the predictions based on bosonization and find reasonable agreement. These results confirm the 
validity of the Tomonaga-Luttinger liquid superlattice description of these systems. 

PACS numbers: 71.10.Pm, 75.45.+J, 75.60.Ej 
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I. INTRODUCTION 

The synthesis of new materials which can be described 
in terms of spin chains or spin ladders has revived the 
study of quantum spin systems in one dimension in the 
last few years i These systems are known to have sur- 
prising features such as the presence of a gap in isotropic 
Heisenberg chains of integer spins and its absence when 
the spins are half-integersiSiSii^Si 7 . More recently, many 
experimental and theoretical results have shown the pres- 
ence of plateaus, in which the magnetization in an exter- 
nal magnetic field is quantized to fractions of the satu- 
rated value &2ii2iiLi2ii^ These systems have special spa- 
tial structures, such as p-merization or ladder geometry, 
which are responsible for the appearance of the magneti- 
zation plateaus. Oshikawa, Yamanaka and Affleck^ de- 
rived the condition p(S — m z ) = integer, necessary for 
the appearance of the magnetization plateaus in ID sys- 
tems. Here, p is the number of sites in the unit cell of 
the magnetic ground state, S is the magnitude of the 
spin and m z is the magnetization per site (taken to be 
in the z-direction). The plateau state can be viewed as 
a gapped state with nonzero magnetization, the integer 
spin chains being a special case where p = 1 and m z = 0. 

Other types of spatial structures, such as quasi- 
periodic couplings, 15 an inhomogeneous magnetic field 16 
or a superlattice structured can also give rise to mag- 
netization plateaus. The superlattice case with periodic 
boundary conditions was studied by two of us in a pre- 
vious work. 17 There, we considered a spin superlattice 
(SS) composed of a repeated pattern of two long and dif- 
ferent spin-i XXZ chains. This model can be viewed as 
the limit of p-merized chains when the number of sites 
per unit cell is very large. Magnetization plateaus were 
found, with magnetization values that depend on the 
relative sizes of the sub-chains, in accordance with the 
condition of Ref. The determination of the width 
and magnetization values of the plateaus relied on the 
Bethe Ansatz exact solution of the XXZ chain^ The 
low-energy properties in the gapless regions away from 
the plateaus, however, could be described by bosoniza- 



tion in terms of a Tomonaga-Luttinger liquid superlattice 
(TLLS)-iS This type of system is of potential interest in 
nanoelectronic applications, where nanowire superlattice 
structures have been built with semiconducting carbon 
nanotubest 2 ^ A TLLS would be obtained if the metallic 
analogue could be manufactured. 

In this work, it is our purpose to check numerically 
the analytical predictions based on the Bethe Ansatz 
and bosonization using the density matrix renormaliza- 
tion group (DMRG)p**2i In particular, we calculate the 
magnetization curve, characterize its plateaus and deter- 
mine the effective Tomonaga-Luttinger liquid parameters 
for spin superlattices. The unit cell of each SS consists 
of two 5=1/2 XXZ sub-chains with different anisotropy 
parameters Aa and sizes L\ (A = 1,2) (but the same 
planar coupling), as shown schematically in Fig.^ The 
Hamiltonian of the SS is written as a sum over N c unit 
cell Hamiltonians 



N c -l 

H ss = H c (j,A u L 1 ,A 2 ,L 2 ) 

3=0 



(1) 



Each unit cell, on the other hand, consists of a total of 
L\ + L 2 — 2 = L c bonds. The first L\ — 1 bonds have 
anisotropy parameter Ai and the following L 2 — 1 bonds 
have anisotropy parameter A2, with Hamiltonian 



Lx-l 



H c (j,A 1 ,L 1 ,A 2 ,L 2 ) = H(jL c + n,A 1 ) 

71=1 
L 2 -l 

+ ]T H {jLc + Li +n- 1,A 2 ), 



where 



H (n, A) — S%S*+i + S^Sf 



n+l 



In the last expression, S%, S% and S z are spin-i operators 
at the n-th site. The total number of sites in the SS is, for 
open boundary conditions, L = (Li + L 2 — 2) 7V C + 1 and, 
for periodic boundary conditions, L = (Li + L 2 — 2) N c . 
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Figure 1: Spin superlattice structure with two unit cells 
(N c = 2), L\ = 4, and L2 = 3. The solid bonds correspond to 
Ai and the dashed ones to A2. The total number of sites is 
L = (Li + L2 — 2) N c + 1 = 11 (open boundary conditions). 



We have set the transverse coupling to 1 to set the energy 
scale. Since we will be interested in increasing both L\ 
and Li while keeping their ratio fixed, we define P as 
the greatest common divisor of L\ and L2 and study the 
behavior of the system as P grows. Note that each sub- 
chain contains L\ — 1 bonds and L\ spins. The spins 
at the boundaries of the sub-chains, however, should be 
viewed as belonging to either sub-chain (see Fig. ^) . Of 
course, in the limit of large L\ in which we are interested, 
this ambiguity at the boundaries is immaterial. Finally, 
the total Hamiltonian has an external magnetic field h 
applied along the anisotropy z-axis 



Hi 



H. 



ss 



L 

£ 

71=1 



hS' 



(2) 



Let us first recall some known results on the XXZ 
model, which is the basic building block of the SS. Using 
the Bethe Ansatz, Yang and YangiS found the exact so- 
lution of the one-dimensional 5 = 1/2 anisotropic (XXZ) 
Heisenberg model. They showed that the model may ex- 
hibit three phases, according to the value of the Ising 
anisotropy A: a ferromagnetic (FM) phase for A < —1, 
a Neel antiferromagnetic (AFM) phase with a spin gap 
for A > 1, and a gapless (critical) phase for —1 < A < 1. 
The low-energy properties of the gapless phase can be de- 
scribed in terms of a Tomonaga-Luttinger liquid^ with 
velocity u and interaction parameter 



II. MAGNETIZATION PLATEAUS 

In the Hamiltonian (J2J, the magnetic field couples to 
a conserved quantity S% ot = J2 n &n- Thus, to obtain the 
magnetization curve we only need the ground state en- 



ergy at h = 0, E{Sl 



0), in each of the subspaces 



with fixed total spin projection 5f ot £ {0, 1, . . . , L/2}. 
Then, we can readily obtain the energy in a finite mag- 
netic field h through the relation E(Sf ot , h) — E(S^ ot ,0) — 
hSf 0t , from which we can construct the magnetization 
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curve i«» 

Since the DMRG is more precise and computationally 
faster with open boundary conditions the magnetization 
curve was calculated this way. We considered lattice 
sizes up to 160 sites keeping up to m = 150 states per 
block. The discarded weight was kept around 10~ 12 . On 
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Figure 2: Magnetization density of a spin superlattice with 
Ai = 0.5, A 2 = 5, Li = 2P, L 2 = 3P {£ = 3/2), and 
N c = 3. The short-dashed line is for P = 3 (L = 40), and 
the long-dashed one for P = 5 (L — 70). The solid line was 
obtained analytically for a spin superlattice with long unit 
cells (P > l)pii 



the other hand, in order to compare with the analytical 
predictions of the Tomonaga-Luttinger liquid parameters 
which were obtained with periodic boundary conditions, 
we have also analyzed SS's with the latter boundary con- 
ditions. In order to obtain a comparable accuracy, we 
considered chains with up to 100 sites with up to m = 600 
states per block. The truncation errors were below 10~ 9 . 

Fig. |21 shows the magnetization density m (magnetiza- 
tion per site normalized to the saturation value) of a SS 
as a function of the external magnetic field for Ai = 0.5 
and A2 = 5. We set N c = 3, since we have observed 
that the magnetization curve is fairly insensitive to the 
number of unit cells. The sub-lattice sizes were chosen 
such that L2/L1 = £ = 3/2. Two lattice sizes are shown 
in the figure: L = 40 and L = 70. They correspond to 
L\ = 2P and L2 = 3P, with P = 3 and P = 5, respec- 
tively. For both sizes the magnetization density has a 
plateau at m = 0.4. 

The continuous line in Fig. [21 corresponds to the case 
where we consider two long sub-lattices with a fixed size 
ratio (P> 1). In this case, the magnetization density of 
the SS is given bjsii 



L1TO1 + L 2 W2 



(3) 



where m\ is the magnetization per site of sub-chain A. 
The magnetization density of each sub-lattice was ob- 
tained from the Bethe Ansatz solution in Ref. UtI 

It can be seen that the three curves have a magnetiza- 
tion plateau at m = 0.4. This has been shown to corre- 
spond to mi = 1 (FM phase) and m-i — (AFM phase) in 
Eq. ©. 17 The plateau is a result of the "spin incompress- 
ibility" of both sub-chains: sub-chain 1 is magnetized at 
saturation, whereas sub-chain 2 has a spin gap. However, 
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3.90474 


2.55006 
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1.44854 


3.51944 


2.07089 


100 


1.47166 


3.36001 


1.88835 


130 


1.48329 


3.27965 


1.79635 


160 


1.48900 


3.23373 


1.74473 


extrapolated 


1.495 


3.146 


1.667 


prediction 


1.5 


3.121 


1.621 



Table I: The saturation field h s , the critical field h c and the 
width of plateau Tnt at m = 0.4 as a function of lattice size. 
For all sizes Ai = 0.5, A 2 = 5, Li = 2P, L 2 = 3P (£ = 3/2), 
and N c — 3. The second to last line is the VBS extrapolation 
and the last one is the prediction of Ref. UtI 



the magnetic field width of the plateau is dependent on 
the sub-chain sizes L\ and L 2 . In fact, the fields at each 
end of the plateau are approached asymmetrically as the 
system grows. The critical field h c (right-hand side of 
the plateau) shows a larger finite size error than the sat- 
uration field h s (left-hand side of the plateau). 

The behavior of the saturation field h s , critical field h c 
and the width of the plateau T jvt = h c — h s as a func- 
tion of the lattice size L is shown in Tabled Again we 
focus on the plateau at m = 0.4 and £ — 3/2 (L\ = 2P, 
L 2 = 3P, and N c = 3), Ai = 0.5, and A 2 = 5. We 
have found the thermodynamic limit of these quanti- 
ties through numerical extrapolation using the Vanden 
Broeck-Schwartz (VBS) algorithm^ which is shown in 
the bottom line. The saturation field increases with the 
lattice size, but slowly at large L, as seen in the first col- 
umn of Table [U Extrapolation to infinite L yields the 
value h s — 1.495. As we have seen in connection with 
Eq. @, we have toi = 1, i.e., the sub-lattice 1 is totally 
magnetized. As shown in Ref. 0, the saturation field is 
the field at which the corresponding homogeneous chain 
with the same anisotropy parameter reaches saturation. 
For Ai = 0.5, this happens at h = 1 + Ai = 1.5, very 
close to the numerical value of h s and compatible with 
our interpretation. 

In the second column of Tabled it can be seen that the 
critical field h c decreases with the lattice size. The con- 
vergence is slower than for h s due to the smaller deriva- 
tive of the magnetization as a function of the field as the 
plateau is approached from above, as is apparent in the 
analytical result of Fig. [21 The extrapolated infinite size 
value of the critical field is h c = 3.146. This should be 
compared with the gap of a homogeneous lattice with 
anisotropy parameter A 2 = 5, which is h = 3.121. This 
value is close to the numerically determined critical field. 
This agreement is again consistent with the sub-lattice 
2 with anisotropy parameter A 2 = 5 being in an AFM 
spin-gapped phase. The width of plateau Tnt = h c — h s 
as a function of L is shown in the third column of Table[U 
The extrapolated infinite size limit of the plateau width 
is = 1.667. 
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Figure 3: The magnetization profile of a spin superlattice at 
the magnetization plateau of m = 0.4. The parameters are 
Ai = 0.5, A 2 = 5, Li = 14, and L 2 = 21 (i = 3/2, L = 100). 
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Figure 4: The width of nontrivial plateau Tnt as a function 
Ai. In this case A 2 = 5, L\ = 6, L 2 = 15 (i = 5/2, m = 
(1 + ey 1 = 2/7), N c = 3, and L = 58. 



In Fig. [3J the magnetization profile at the plateau of 
to = 0.4 is shown. The parameters used are L = 100 
(Li = 2P, L 2 = 3P, P = 7, and N c = 3), I = 3/2, 
Ai = 0.5, and A 2 = 5. We can see that, indeed, in 
the sub-lattices with anisotropy parameter Ai = 0.5, the 
spins are fully polarized, whereas the sub-lattices with 
anisotropy parameter A 2 — 5 are antiferromagnetically 
ordered. 

The overall picture resulting from Figs. [21 and [3 and 
Table[I]is thus compatible with the SS having a nontrivial 
plateau at to = 1/(1 + £), in which one sub-chain is fully 
saturated with toi = 1, whereas the other is in a spin- 
gapped AFM phase with to 2 = 0iL2& 

The width of the plateau of a SS with A 2 = 5 is shown 
in Fig. 21 as a function of the anisotropy parameter Ai. 
The other parameters are L\ = 2P, L 2 = 5P (£ = 5/2, 
to = (1 + fT 1 = 2/7), P = 3, 7V C = 3 5 an d L = 58. We 
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Hamiltonian of TLLS is given by 
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Figure 5: The width of the plateau at zero magnetization as 
a function of the lattice size for a SS with Ai = 2, A2 = 5 and 
e = 3/2(Li = 2P, L 2 = 3P, and N c = 3) (closed triangles). 
For comparison, we also show the size dependence of the spin 
gap of a homogeneous chain with anisotropy parameter A = 2 
(open diamonds). The extrapolated infinite size limits of the 
two curves (closed circles) coincide within the numerical error. 



observe that the width of the plateau decreases linearly 
with Ai and vanishes at Ai = 2.741. The magnetization 
curve of this SS does not have plateaus for Ai > 2.741. 
This linear dependence is expected. As we have seen, the 
plateau width is given by Tnt — hf(A 2 = 5) — 1 — Ai, 
where is the gap of the corresponding homogeneous 
spin chain in the AFM phased The slope of the line in 
the Fig. ^ is —0.966, close to the expected value of — 1. 

The magnetization curve in Fig. |3 has one nontrivial 
plateau at m = 0.4 for the parameters Ai = 0.5 and 
A 2 = 5. When both Ai > 1 and A 2 > 1, however, 
a new plateau at m — emerges. In this case, both 
sub-chains 1 and 2 are in a spin-gapped AFM phase. 17 
For a SS with parameters Ai = 2, A 2 = 5, and i — 
3/2 (£1 = 2P, L 2 = 3P, and N c = 3), the width of 
the plateau at m = is shown in Fig. [S] as a function 
of 1/L. Using the VBS algorithm, we determined the 
infinite size limit of the plateau width to be Ft — 0.198. 
According to the analysis of Ref. 0, this value should 
be given by the gap of a homogeneous spin chain with 
anisotropy parameter Ai = 2. For comparison, we also 
show in Fig. |S] the value of the latter gap as a function 
of system size. The extrapolated gap size is Tt = 0.196, 
while the value obtained from the exact solution is Tt = 
0.19842. Indeed, the extrapolated values coincide within 
the numerical error. 



III. GAPLESS REGION: THE EFFECTIVE 
TOMONAGA-LUTTINGER PARAMETERS 

The low-energy properties of the SS away from the 
plateaus can be described in terms of a TLLS^* The 



Hllsl 



1 

2^ 



dx \ u{x)k(x) (d x ey 



u(x) 



K(x) 

where d x Q is the momentum field conjugate to $: 
[&(x),dyQ(y)} — iS(x — y). The fields 4> and 6 can be 
related to the spin density operators^ In the Hamil- 
tonian l|4jl. we have introduced the sub-chain-dependent 
parameters u(x) and K(x). For x in the sub-chain A, one 
has K(x) = K(J, Ax, h) and u{x) — u(J, A\, h), i.e., the 
usual uniform Tomonaga-Luttinger parameters for each 
sub-chain, which can be obtained directly from the Bethe 
Ansatz solutionii2*S Of course, this effective low-energy 
description is valid asymptotically in the limit of very 
long sub-chains. Using periodic boundary conditions and 
diagonalizing the Hamiltonian we find that the low 
energy properties of the SS are determined by just a few 
effective parameters ^LA£ These parameters are the effec- 
tive velocity c and the effective correlation exponents K* 
and K, which are given byi£*i& 



\Jl + T]lu\ju2 + (£ui/u2) 2 



(5) 



K* 



1 + rj£ui/u2 + {lux/u-z) 
_i L 1 Ul 



K = f(l/K u l/K 2 ), 



f(K 1: K 2 ), (6) 



(7) 



where 77 = K 1 /K 2 + K 2 /Ki. Clearly, {c,K*^K) -> 
{u 2 , K 2 ,l/K 2 ) as I -» 00, and (c,K*,K) -> 
(ui,Ki,l/K{) as I — > 0, as expected. The important 
feature to notice in Eqs. JSJ, iJBJ, and J7J is the fact that 
the effective SS parameters represent a certain weighted 
average of the individual sub-chain velocities and corre- 
lation exponents. This weighted average is induced by 
the superlattice structure and is a feature ubiquitous in 
TLLS's-iLia 

It is straightforward to extract from the Hamiltonian 
Q the finite-size spin gap of the system. It is given by 



E{Sf ot = l,h = 0)-E(Sf ot =0,h = 0) = 



2K*L 



(8) 



Thus, from the scaling of the spin gap with the system 
size, we can verify the predictions of Eqs. JSJ and l@ for 
the effective Tomonaga-Luttinger parameters. 

Fig. shows the numerically determined ratio c/K* 
for a SS with Ai = 0.4 and A 2 = 0.8 as a function of 
t We used N c = 4 and L = 100, 76, 100, 52, 76, 100, for 
t = 1/2, 3/4, 5/4, 3/2, 5/2, 7/2, respectively. For compar- 
ison, we also show the TLLS prediction obtained from 
the ratio of Eqs. JSJ and © and from the known val- 
ues of u\ and K\ for homogeneous chains. We can see 
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Figure 6: c/K* as a function I for a SS with Ai = 0.4 
and A2 = 0.8. The continuous line was obtained using 
bosonization»'£ The closed dots correspond to the DMRG re- 
sults. 

that there is reasonable agreement, with slightly larger 
discrepancies at larger I. The ratio u/K for homoge- 
neous chains with anisotropy parameters Ai = 0.4 and 
A 2 = 0.8 are equal to U\/K\ — 1.57 and U2/K2 — 2.33, 
respectively. c/K* interpolates smoothly between u\ /K 1 
and U2I K2 as £ increases, a manifestation of the spatial 
averaging due to the superlattice structure. We believe 
the small discrepancies between the curves in Fig. [6] are 
due to the finite sizes of the sub-chains. We recall that 
the TLLS predictions are expected to hold asymptoti- 
cally for very long sub-chains. For a gapless phase, the 
inhomogeneities created by the boundaries between sub- 
chains will give rise to Friedel oscillations which die out 
only as power lawsi^i These disturbances are expected 
to give rise to finite-size corrections to the TLLS pre- 
dictions. We stress, however, that although the TLLS 
analysis predicts a sort of weighted average for the de- 



pendence of c/K* on £, the detailed form of this average 
is highly nontrivial. Yet, precisely this non-linear depen- 
dence is strikingly confirmed by the numerical data. We 
consider this as a stringent test of the predictions of the 
theory. 



IV. CONCLUSIONS 

In summary, we have used the finite size DMRG 
method with open and periodic boundary conditions 
to study spin superlattices made up of a periodic ar- 
rangement of two XXZ chains with different parame- 
ters and sizes. We confirmed previous analytical predic- 
tions of a nontrivial plateau in the magnetization curve 
at m = 1/(1 + £), where £ is the relative size of the 
sub-chains. When both anisotropies are larger than 1, 
we have also confirmed the expected trivial plateau at 
m = 0. The nontrivial plateau width was shown to ap- 
proach the asymptotic value of a superlattice with long 
sub-chains. Moreover, the magnetization profile was seen 
to be in accord with the analytical predictions of one 
sub-chain being saturated and the other one being AFM 
ordered. Finally, we found fairly good agreement in the 
gapless region with the results of a Tomonaga-Luttinger 
liquid theory as applied to a superlattice structure. This 
agreement confirms the non-trivial prediction of this the- 
ory for the way the individual sub-chain properties are 
averaged over in the effective low-energy description of 
the superlattice. 
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